Average Patterns of Spatiotemporal Chaos: A Boundary Effect 
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Chaotic pattern dynamics in many experimental systems show structured time averages. We suggest that 
simple universal boundary effects underly this phenomenon and exemplify them with the Kuramoto-Sivashinsky 
equation in a finite domain. As in the experiments, averaged patterns in the equation recover global symmetries 
locally broken in the chaotic field. Plateus in the average pattern wavenumber as a function of the system size 
are observed and studied and the different behaviors at the central and boundary regions are discussed. Finally, 
the structure strenght of average patterns is investigated as a function of system size. 
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Experimental studies of the chaotic pattern dynamics in 
Faraday waves in rotating thermal convection and 

in electroconvection [Q], reveal that spatio-temporal complex 
patterns can have surprisingly ordered time averages. The 
form of these average patterns (square, circular, hexagonal) 
is determined by the underlying symmetry imposed by 
the boundary conditions. Although the instantaneous patterns 
fluctuate chaotically, they are biased towards the average pat- 
tern because they have short-lived patches spatially in phase 
with this average. This phase rigidity seems to come from 
the boundaries, and quantization effects appear due to the fi- 
nite size of the container The amplitude of the time-averaged 
pattern depends on the system size and control parameters. 
It is strongest near the sidewalls, and decays with increasing 
distance from the sidewalls and with increasing fluctuations 
about the ordered averaged state. For very large containers 
the ordered average pattern exists only near the sidewalls. 

Given the general features of average patterns suggested by 
experiments, it seems surprising that their possible existence 
and characterization have not been addressed within the stan- 
dard model equations displaying spatiotemporal chaotic states 
10] . One possible reason for this is that periodic boundary con- 
ditions are usually considered in theoretical studies. In such 
situation spatial translational invariance homogenizes out any 
time average (unless some unexpected ergodicity breaking 
takes place). Boundary conditions breaking translational sym- 
metry, as in the experiments, are thus needed to obtain non- 
trivial average patterns. Motivated by this fact, we here con- 
sider the Kuramoto-Sivashinsky equation, one of the proto- 
type equations showing spatiotemporal chaos, in bounded one 
and two dimensional domains. We show that ordered aver- 
age patterns do appear, despite the strong fluctuations, and we 
discuss the universal aspects of wave number selection and 
amplitude variations. More directly, our analysis of average 
patterns may be relevant and suggestive for experiments on 
phase turbulence in convection cells fluids flowing down 
an inclined wall [^], and flame front propagation [ 10,11]]. 

The Kuramoto-Sivashinsky equation [pljl^] is perhaps 
the simplest partial differential equation exhibiting spatio- 
temporal chaos. The equation in one dimension has the form: 



hf = -K 



(1) 



where h ~ h{x, t) is a real function, x E [0, L], and the sub- 
scripts stand for derivatives. In two dimensions the spatial 
derivative is replaced by a gradient and the second derivative 
by a Laplacian. The only control parameter for the equation is 
the length of the domain L; prefactors to the terms in Eq. (|l]) 
can be scaled out. An equivalent equation for u = hx can be 
obtained by taking the derivative of Eq. dl]) with respect to x. 
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Eq. (jlj) possesses translational symmetries (h h + ho, x 
X + xq) a reflexion symmetry (h ^ h, x ^ ~x), and an 
infinitesimal Galilean symmetry (x — > a; + 2vt, h ^ h + vx). 
Eq. (H) is also invariant under translations (u u, x x + 
xq) and under a Galilean symmetry (x — + x+2vt, u — > u+v). 
A different reflexion symmetry is valid in this case (u ~* ~u, 
x —x). 

The stability of the laminar solution /i = (u = 0) is an- 
alyzed by linearizing Eq. (|l]) [Eq. (^]. For commonly used 
boundary conditions the growth rate /i for the Fourier mode 
of wavenumber fc is /i = fc^ — /c'*. In two dimensions fc^ 
is replaced by |kp. The laminar solution is unstable for all 
modes within < A; < 1. The fastest growing mode has 
a wavenumber fcg = 1/ \/2 corresponding to a wavelength 
Ao = 2\/2vr « 8.9. The wavelength Aq serves as a basic 
length scale, and the system size L is naturally measured in 
units of this scale, L/Xq, which is called the aspect ratio. Be- 
yond the linear range, the nonlinear term becomes important 
and produces growth (linear in time) of the mean value of h, 
while the mean value of u saturates. For L large enough to 
permit a sufficient number of unstable Fourier modes, the so- 
lution exhibits spatio-temporal chaotic behavior that can be 
associated with a disordered evolution of a cellular pattern. 

Many studies have been devoted to the bulk behavior of the 
Kuramoto-Sivashinsky system In relation to average pat- 
terns however the boundaries are of paramount importance, 
as discussed above. Here, we consider two types of bound- 
ary conditions. One of them is the rigid boundary conditions, 
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where 

u(0, t) = u{L, t) = u^{Q, t) = u^{L, i) = , (3) 

or equivalently, 

t) = h,{L, t) = h,,{0, t) = h,,{L, t)=0. (4) 

Our other choice of boundary conditions is 

u{0, t) = u{L, t) = u,,(0, t) = <) = , (5) 

or equivalently, 

h,{0, t) = K{L, t) = /i,,,(0, t) = /i,,,(L, t) = , (6) 

which we call stress-free boundary conditions, with reference 
to similar conditions in hydrodynamics. 

We integrate the Kuramoto-Sivashinsky equation using ex- 
plicit finite-differences of first order in time, second order in 
space for the linear terms, and fourth order in space for the 
nonlinear term. The time step is chosen sufficiently small to 
avoid any spurious behavior The number of grid points used 
is 128 in one dimensional simulations and 64 x 64 in two 
dimensions. In all cases, the simulations were started from 
random initial conditions. 
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FIG. 1. Results from simulations of the one-dimensional Ku- 
ramoto-Sivashinsky equation with stress-free boundary conditions. 
System size L = 70. a) A characteristic front (p at time t — 1000. b) 
The time average of the front ip (continuous line) and of u (dashed). 

We are interested in the average pattern of u and the average 
pattern of the front cp of h, 

1 

ip{x,t) = h{x,t) — - / h(x,t)dx . (7) 
L Jo 



To optimize the measurements of the average, the sampling 
was first started well beyond the initial transient behavior For 
the system sizes considered, the typical transient time was 
limited to approximately 20 time units, and we discarded the 
first 100 time units. Then, averages where taken from config- 
urations sampled every 5 time units. A total of 10, 000 con- 
figurations per run were included, and further average over 
10 runs with independent random initial conditions was per- 
formed. This is a large sample, but was necessary to compen- 
sate for the slow convergence of the averages produced by the 
long-range time correlations present in the KS equation. 
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FIG. 2. The same as Fig. |^ but for rigid boundary conditions 

in both one and two dimensions and with both rigid and 
stress-free boundary conditions we obtain non-trivial and 
ordered time-averaged patterns from the spatio-temporally 
chaotic evolution (Figs. |l}§), emphasizing that the forma- 
tion of average patterns in spatio-temporal complex systems 
is general despite the presence of very large fluctuations. The 
presence of boundaries breaks the translational symmetry of 
the equations. The boundary conditions (|^-(|1) respect how- 
ever its reflection symmetries (for reflexions with respect to 
the center of the domain). As in the experiments, here we find 
that the average patterns display these remaining symmetries. 
In the two-dimensional case the average pattern recovers also 
the square symmetry of the integration domain (Fig. ||). Ex- 
cept for the one-dimensional case with stress-free boundary 
conditions (Fig. |l]), an overall parabolic profile of the average 
front is obtained (see Figs. For the derivative u a mean 
slope is obtained [|r3]]. This parabolic profile is a peculiarity 
of the Kuramoto-Sivashinsky equation. It can be removed by 
considering the second derivative of the front ip instead of ip 
itself; for this variable (and for the Laplacian in two dimen- 
sions) the discussion for all the cases is very similar to the 
one-dimensional stress-free situation, that we address in fur- 
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ther detail in the remaining of the paper. 

Fig. ^ shows the average patterns for L — 60 and L = 100. 
The number of oscillations increases with the size of the sys- 
tem, although only those close to the boundaries are large. 
Furthermore, the distance between consecutive maxima is 
close (but not equal, see below) to the characteristic length 
scale Ao ~ 8.9. Similar observations were done in the ex- 
periments referred to at the beginning. In Fig. ^ the number 
of local maxima N in the average front is shown for increas- 
ing values of the aspect ratio L/Xq. N is written in terms 
of the average distance A between two consecutive minima, 

= L/\. The line where A = Aq is also indicated in Fig. ^ 
Plateaus are obtained at every integer N between 6 and 11 for 
system sizes L between 59 and 112. The average distance A is 
consistently larger than Aq. Consistent deviations (positive or 
negative) are also known from the Faraday wave experiments 




FIG. 3. Results from simulations of the two-dimensional Ku- 
ramoto-Sivashinsky equation with stress-free boundary conditions. 
System size L x L = 70 x 70. Left: A characteristic instantaneous 
image of the front tfi{t = 1000), different values of are coded with 
different gray levels. Right: The average of V^^p, showing the square 
symmetry. The average of the Laplacian instead of the front itself is 
shown to eliminate the dominant parabolic shape of the average, thus 
improving the visibility of the sidewall oscillations. 

If only the central region is considered, the plateaus fall off. 
More specifically, consider the 'central region', defined as the 
domain ranging from the second local minimum to the sec- 
ond last local minimum of </? (see Fig. The rest of the the 
pattern is thus considered the 'boundary region' . We now de- 
termine the average distance Ac between consecutive minima 
in the central region for various system sizes L, and find the 
number of maxima L/\c characteristic for the central region. 
The results are shown in Fig. |[ Intriguingly, the plateaus now 
fall off, an observation also done in experimental studies of 
the central region [|l]]. In order to explain this 'fall off effect, 
we determine the average distance Af, between minima in the 
boundary regions. Over the entire range of system sizes con- 
sidered, these distances changes very little, not more than 4%, 
so that to a first approximation we can consider Af, indepen- 
dent of L. For the central region we now have 



A: = (^-4)7r- 
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(8) 



plateau characterized by N. Thus an almost constant value of 
Af, serves as a generic explanation for the generally observed 
fall off of the plateaus. The overall picture is that when L 
is increased the total number of oscillations tends to remain 
constant, as well as Af,, so that Ac increases. This situation 
continues until the local wavelength in the central region is 
far enough from Ao, moment at which a new oscillation is ac- 
commodated and a jump in N occurs. 




The last approximation is valid for Af,/L small. For Af, con- 
stant, it is seen that L/Xc falls off as ~ within a given 



FIG. 4. a) A central horizontal cut of the configuration shown in 
the left panel of Fig. ^ b) Solid: a central horizontal cut of the time 
average of the front tp for the same parameter values; dashed: the 
same cut for the time average of ipx ■ 

From Fig. ||it is clear that the amplitude A{x) of the aver- 
age pattern in general decays with increasing distance from 
the boundaries |]l4|]. Experimental studies show the same 
behavior [|l]-0]. To quantify this observation, we consider 
the spatial average (^^) = A^dx. The variation of 

A-rms = VT^^ system size is shown in Fig. ^ show- 
ing a power-law dependence as Arms ^ L^^l"^ . We explain 
this fact by noting that Fig. ^indicates that A^ms receives its 
largest contribution from the boundaries, so that the integral 
in the definition of (^^) becomes a constant for system sizes 
larger than the boundary region. Thus the factor L^^ in the 
definition of (^^) becomes the dominant L-dependence thus 
providing the observed behavior of Arms- 

In conclusion, we have established the formation of ordered 
time-averaged patterns in the Kuramoto-Sivashinsky equa- 
tion, in one and two dimensions, and with rigid as well as 
stress-free boundary conditions. The average pattern recov- 
ers the symmetries which are respected by both the equation 
and the boundary conditions. The amplitude is strongest at 
the boundaries and decays with increasing distance to them. 
The law of decay has been found and explained. We have 
determined the selected wavelength A, its variation with sys- 
tem size L, and interpreted the different behavior between 
the central and boundary regions. Most of these observa- 
tions are also found in experimental systems for which the 
Kuramoto-Sivashinsky equation does not apply, thus indicat- 
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ing its generic, mainly geometrical, origin: What is relevant 
for these phenomena to occur is the occurrence of strong 
enough chaotic fluctuations in the presence of non-trivial 
boundaries. 
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FIG. 6. -I-: The number of maxima A'^ = L/X 'm the average front 
versus the aspect ratio L/Aq, A being the average distance between 
consecutive minima in the entire region, o: the number of maxima 
in the central region, given by L/Ac. The solid line corresponds to 
A = Ao. 
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FIG. 5. Average patterns of the front <^ for the one-dimensional 
Kuramoto-Sivashinsky equation with stress-free boundary condi- 
tions. &) L = 60. b) i = 100. The arrows indicate the 'central 
region' , as defined in the text. 



FIG. 7. Variation of the root mean square amplitude Arms over 
the range of system sizes 59 < L < 107. The solid curve is the 
function C /\/L, with C fitted to the data 
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